Kappa-deformed random-matrix theory based on 
Kaniadakis statistics 



A.Y. Abul-Magd 1,2) and M. Abdel-Megeed 2) 
V Faculty of En gineering, Sinai University, El-Arish, Egypt 
2 ) Faculty of Science, Zagazig University, Zagazig, Egypt 

January 4, 2012 



Abstract 

We present a possible extension of the random-matrix theory, which 
is widely used to describe spectral fluctuations of chaotic systems. By 
considering the Kaniadakis non-Gaussian statistics, characterized by the 
index k (Boltzmann-Gibbs entropy is recovered in the limit k — > 0), we 
propose the non-Gaussian deformations (k 7^ 0) of the conventional or- 
thogonal and unitary ensembles of random matrices. The joint eigenvalue 
distributions for the K-deformed ensembles are derived by applying the 
principle maximum entropy to Kaniadakis entropy. The resulting dis- 
tribution functions are base independent as they depend on the matrix 
elements in a trace form. Using these expressions, we introduce a new 
generalized form of the Wigner surmise valid for nearly-chaotic mixed 
systems, where a basis-independent description is still expected to hold. 
We motivate the necessity of such generalization by the need to describe 
the transition of the spacing distribution from chaos to order, at least in 
the initial stage. We show several examples about the use of the general- 
ized Wigner surmise to the analysis of the results of a number of previous 
experiments and numerical experiments. Our results suggest the entropic 
index k as a measure for deviation from the state of chaos. We also 
introduce a K-deformed Porter-Thomas distribution of transition intensi- 
ties, which fits the experimental data for mixed systems better than the 
commonly- used gamma-distribution. 



1 Introduction 

Quantum systems whose classical counterparts have chaotic dynamics are often 
modelled in terms of random- matrix theory (RMT) [TJ [2J. This is the sta- 
tistical theory of random matrices H whose entries fluctuate as independent 
Gaussian random numbers. Dyson [3] showed that there are three generic en- 
sembles of random matrices, defined in terms of the symmetry properties of 
the Hamiltonian. The most popular of these is the Gaussian orthogonal ensem- 
ble (GOE) which successfully describes many time-reversal-invariant quantum 
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systems whose classical counterparts have chaotic dynamics. Systems without 
time- reversal symmetry are modelled by the Gaussian unitary ensemble (GUE). 
Balian 4 derived the weight functions for the random-matrix ensembles from 
Jaynes' principle of maximum entropy (MaxEnt) [5]. He applied the conven- 
tional Shannon entropy to ensembles of random matrices with distribution P(H) 
and maximized it under the constraints of normalization of P{H) and fixed mean 
value of Tr[WHj, where is the Hermitian conjugate of H. The latter con- 
straint ensures basis independence, which is a property of the trace of a matrix. 
The resulting expression 



is supposed to account for all possible level statistics of various chaotic systems. 

For most systems, however, the phase space is partitioned into regular and 
chaotic domains. These systems are known as mixed systems. Attempts to 
generalize RMT to describe such mixed systems are numerous; for a review 
please see [2J. MaxEnt can contribute to the generalization of RMT by either 
introducing additional constraints or modifying the entropy. An example of the 
first approach is the work of Hussein and Pato [TJ, who use MaxEnt to con- 
struct "deformed" random-matrix ensembles by imposing different constraints 
for the diagonal and off-diagonal elements. The second approach to include 
mixed systems in RMT is to apply MaxEnt to entropic measures other than 
Shannon's. Indeed, there is no systematic way of deriving the "right" entropy 
for a given dynamical system. Different entropic forms have been introduced 
instead of Shannon's entropy, which would generalize the classical statistical 
mechanics of Boltzmann and Gibbs. During the last two decades a great inter- 
est has arisen in the non-extensive Tsallis statistical mechanics [8], and, more 
recently, on the Kaniadakis extensive generalized power-law statistics [9] moti- 
vated by various restrictions to the applicability of classical, extensive statistical 
mechanics. This has been accomplished with the introduction of an appropri- 
ate deformed logarithm [?]. Recently, the Tsallis formalism has been applied 
to include systems with mixed regular-chaotic dynamics in a nonextensive gen- 
eralization of RMT PHI [TTJ [H EE HI EU- The matrix-element distribution 
P(H) is obtained by extremizing Tsallis' g-entropy rather than Shannon's, but 
again subject to the same constraints of normalization and constant expectation 
value of Tr (WH\. This distribution function, as well as the distributions to be 
obtained for other generalized entropies, depends on the matrix elements in the 
combination Tr (if ^if). Therefore, the random-matrix ensembles described by 
these distribution functions are base independent. 

Another nonextensive approach to RMT [IB] replaces Shannon's entropy in 
the conventional theory by Kaniadakis' K-entropy [9l [TT1 [18l [19] . The Kaniadakis 
entropy is maximized again under the constraints of normalization of P(H) 
and fixed mean value of Tr (H^H). Clearly, the resulting distribution function 
depends on the matrix elements in the combination Tr(_ff 'ffl as in the case 
with Tsallis' entropy. Therefore, the emerging random-matrix ensemble is also 
base independent. 




(1) 
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In this paper, by following our previous results [TB], we adopt the Kaniadakis 
non-Gaussian statistics and study the effects of the K-generalization for RMT, 
especially what concerns its effects on NNS of the eigenvalues of GOE and GUE 
of random matrices. We also investigate the influence of these non-gaussian 
effects on transitions in the mixed system out of the chaotic phase. 

This paper is organized as follows. A brief review of Kaniadakis statistics 
is presented in Sec. 2. In Sec. 3, we consider the connection between the 
K-generalization of RMT and the transition of a quantum system out of the 
chaotic phase. Section 4 presents a K-generalization of the Wigner surmise 
which is a good approximation for the spacing distribution of chaotic systems. 
A comparison between the predictions of the K-deformed Wigner surmise with 
the results of a numerical experiment is given in Sect. 5. Section 6 presents 
a formula for the transition intensities which is shown to be more consistant 
with experimental data than the gamma distribution which is commonly used 
in data analysis. We summarize our main conclusions in Sect. 7. 



2 Kaniadakis' entropy 

In this section, we present the basic elements of Kaniadakis' K-deformed frame- 
work introduced in [5]. The Kaniadakis non-gaussian statistics in turn is char- 
acterized by the K-entropy that emerges naturally in the framework of the so- 
called kinetic interaction principle, which underlies the non-linear kinetics in 
particle systems. This principle fixes the expression of the Fokker-Planck equa- 
tion describing the kinetic evolution of the system and imposes the form the 
generalized entropy associated with the system. The structure of the ensuing 
K-deformed statistical mechanics has striking similarity with that of special rel- 
ativity, suggesting that it might be relevant for the self-consistent formulation 
of relativistic statistical theory [18l [19] and relativistic Boltzmann kinematics 
[20] with application to, e.g., relativistic distribution of fluxes of cosmic rays 
particle systems [JJ5] , the formation of a quark-gluon plasma [2T] , the relativis- 
tic nuclear equation of state [22] , the relativistic gas in an electromagnetic field 
[2"3"] , and the relaxation of relativistic plasmas under the effect of wave-particle 
interactions [53]. The ft-deformed statistics have also been successfully applied 
in several branches of physics. Examples are the nonlinear kinetics [20], and 
the H-theorem from a generalization of the chaos molecular hypothesis [25] , the 
sensitivity to initial conditions and the entropy production in the logistic map 
[2"S] , personal income distribution for Germany, Italy and UK [27J , the distribu- 
tion of stellar rotational velocity for stars belonging to various classes [2E] and 
systems of interacting atoms and photons [29] . 

From the mathematical point of view, Kaniadakis' statistics is based on the 
K-deformed exponential I9i , which is defined by 



exp K (x) = 1 + k 2 x 2 + xj = exp f — arcsinh nxj . (2) 
The properties of the function exp K (x) have been considered extensively in the 
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literature. We recall briefly that in the n — > limit the function exp K (x) reduces 
to the ordinary exponential, i.e. has the properties exp (x) = exp(a;). For 
x — > 0, it behaves very similarly with the ordinary exponential independently 
on the value of k. It is remarkable that the first three terms of the Taylor 
expansion 

2 3 

ex.p K (x) = l + x + — + (1 - K 2 ) + ... (3) 

are the same as in the ordinary exponential. On the other hand, the most 
interesting property of exp K (x) for the applications in statistics is the power- 
law asymptotic behavior 

exp K (z) ~ |2 ra | ±1/|K| . (4) 

x— > ±oo 

The corresponding logarithmic function ln K (x) is defined as the inverse function 
of exp K (x), namely ln K (exp K (x)) = exp K (ln K (x)) = x, and is given by 

ln K (x) = (5) 
Zk 

Kaniadakis' entropy associated with the K-statistics is obtained by replacing 
the logarithm in the expression for the standard Shannon's entropy by the k- 
logarithm. For a distribution function /, the /t-entropy is given by 

S K = -[dTfln K f (6) 



where dT is the phase-space volume element. Maximizing the K-entropy with 
condition of fixed mean value of a Hamiltonian h yields 

P(H) = lexp K (-/3/0, (7) 

Zk 

where /? is a Lagrange multiplier, and z K is a moralization factor. We see 
that the Gaussian probability curve is replaced by the characteristic power-law 
behavior of Kaniadakis framework and, as expected, the limit k = recovers 
the exponential result. 



3 The ^-deformed RMT 

In this section, we consider a possible generalization of RMT based on an ex- 
tremization of Kaniadakis' K-entropy. We obtain the Kaniadakis entropy for 
the matrix-element probability distribution function using Eqs. © and ©. 
Explicitly, it reads 

S K [k,P k (k,H)) =-± [dH [P k (k,H)] 1+k - [P^H)] 1 -*) 

(8) 
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where k is a parameter with value between and 1. The case of K = 0, corre- 
sponds to the Shannon entropy. Here, a is a real positive parameter. Kaniadakis 
has considered two choices of a, namely a = 1 and a = Z, where Z is a general- 
ized partition function. We here adopt the second choice. The matrix-element 
distribution Pk(k,H) is obtained by extremizing the K-entropy subject to the 
condition of fixed mean value of Tr [H^Hj, One has to obtain the extremum of 
the functional 

F K = S K -r lK J dH P K ( K , H)Tt (H t H) , (9) 
where rj K is a Lagrange multiplier. One arrives to the following distribution 

P k (k, H) = I cxp K [-77 K Tr (H T H)] , (10) 

Z = JdH exp K [-77 K Tr {H T H)] . (11) 

Accordingly, the K-deformed matrix-element distribution Pk(k, H) is obtained 
by replacing the exponential function in Eq. (TT]) by the ^-exponential given by 
identity ©. One arrives to the following distribution 

ifc(K,ff) = i exp K [-r ?K Tr(fft jff )] ! ( i2 ) 

in which 

Z K = J dH cxp K [-77k Tr H)} bv. (13) 
is a normalization factor defined by 

1 \ / \ ™/ 2 F f — 



where 



Z K = 2-^-^1* J <fWp (--J-arcsinh K ry K r 2 ) = 2-W"- 1 )/ 4 ^ 
for k < 2/n. Here we define a K-deformed gamma function by 



(14) 



i \ i r fnl r (— -) 

T K (n) = I X™" 1 exp ( arcsinh kx ) dx = , \ 2k (15) 



K 



(2 K ) 1+n r(i + i + f) 



A similar result is obtained by Kaniadakis in evaluating an analogous integral 
that appears in his treatment of Brownian particles (Eq. (71) of Ref. [9]). 

We now calculate the joint probability density for the eigenvalues of the 
Hamiltonian H. With H = U~ 1 EU, where U is the global unitary group, we in- 
troduce the elements of the diagonal matrix of eigenvalues E — diag(2?i, • • • , En) 
of the eigenvalues and the independent elements of U as new variables. Then 
the volume element dH has the form 

dH =\A N (E)\dEdfj,(U), (16) 
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where Ajy (E) = Yl n>m (E n — E m ) is the Vandermonde determinant and dfi(U) 
the invariant Haar measure of the unitary group [TJ [2] . The probability density 
Pk(k,H) depends on H through TrLH'H} and is therefore invariant under 
arbitrary rotations in the matrix space. Integrating over the " angular variables" 
U yields the joint probability density of eigenvalues in the form 



P k (k;E 1 ,...,E n ) = C H \E n -E m fexp K 



-Vk 



N 
i=l 



(17) 



where C is a normalization constant. 



4 The K-deformed Wigner surmises 

The calculation of NNS distributions starting from the matrix-element proba- 
bility density function ((T|) is a formidable task. In fact, RMT does not provide 
simple analytical expressions for NNS distributions. There are several elaborate 
approaches to evaluate these distributions T| . The author is not aware of appli- 
cation of these results to the analysis of experimentally observed or numerically 
calculated discrete data. 

In many cases, the empirical data are actually compared to the so-called 
Wigner surmises, which correspond to the NNS distributions of ensembles of 
2x2 matrices. For the three Gaussian ensembles, the Wigner surmises are 
given by 

Pp(s)=a^exp(-b p s 2 ), (18) 

where /3 = 1,2, and 4 for GOE, GUE and GSE, respectively. The factors ap 
and bp are obtained from the normalization condition and the requirement of a 
unit mean spacing, 

/>oo />OC 

/ Pp(s)ds = l, / sPp(s)ds = l. (19) 
Jo Jo 

For GOE, at = vr/2 and b x = tt/4 while, for GUE, a 2 = 32/vr and b 2 = 4/vr. 
The Wigner surmises present accurate approximation to the exact results for 
the case of large N. 

We consider the K-deformed generalization of the Wigner surmises hoping 
that they present an accurate approximation for ft-deformed ensembles of large 
N, which has an equal success to that of the standard Wigner surmises. For 
this purpose we follow the arguments given in the work of Carvalho et al. |28j . 
We replace the exponential function in Eq. (1181) by the K-deformed exponential 
in Eq. ©to obtain 

P K (/3, k- s) = a(j3, k) S ? exp K [-b(/3, k)s 2 ] . (20) 

Again the factors a(/3, k) and 6(/3, k) are obtained from the normalization and 
unit-mean-spacing conditions (I19p so that 

,„ s 26 £ * i (/3, K ) 
a(P,K) = j^-^ (21) 



G 



and 



608, k) 



(«+0 



1 « I 2 



(22) 



We can easily sea that a(l,0) = 7r/2 and 6(1,0) = 7r/4 while a(2,0) = 32/tt 
and 6(0, k) = Therefore, as exp (x) = exp(x), the proposed K-deformed 

Wigner surmises Pk(P, k; s) tend to the conventional ones as n — > 0. Using the 
small argument expansion for exp K (x) © , Eq. (j2"0")) yields 



P k (/?,k; s) = a{P, K ) S P [l + 0(s 2 )] 



(23) 



suggesting that the /c-deformed Wigner surmises have the same power-law level 
repulsion as the original ones. On the other hand, the large-argument behavior 
of exp K (x) in Eq. (j4]) implies the following power-law asymptotic behavior of 
Pk(0, k] s) 



Pk(/3, k; s) oc s 1 2/k , for s > 



1 



V / «6(^, k) 



(24) 



The mean square spacing of the distribution (I2U1) . (s 2 ) = J* °° s 2 P/3(s)ds 7 is 
given by 



<* 2 > 



\^2k ' 4 



4«rfi-£-±^r(i + £ + 2 



(25) 



It is reasonable to suggest that the condition of fixed mean value of Tr(wH) 
implies that (s 2 ) is finite. This restricts the range of allowed values of n to the 
following 

< k < — . (26) 



5 NNS distribution for mixed systems 

Bohigas et al. |32) put forward a conjecture (strongly supported by accumulated 
numerical evidence) that the energy levels of a typical time-reversal invariant 
chaotic systems have the same fluctuation statistics as that of GOE. In partic- 
ular, the NNS distribution is given by the Wigner surmise (|18|) with j3 = 1 

P G oe(s) = |scx P (-J s 2 ) (27) 

In the nearly ordered regime, mixing of quantum states belonging to adjacent 
levels can be ignored and the energy levels are uncorrelated. The level-spacing 
distribution function obeys the Poissonian 

fpoisson(s) = exp(-s) (28) 
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Systems with mixed regular-chaotic dynamics have NNS distribution interme- 
diate between the Wigner and Poisson distributions. Several formulae were 
suggested to interpolate between these two distributions [2J. 

The NNS statistics of mixed quantum systems execute a Poisson-to- Wigner 
transition as the underlying classical dynamics monotonically change from being 
completely integrable to completely chaotic. We explore in this section the 
relevance of the K-deformed Wigner surmise ([20)) for mixed systems. We do 
this phenomenologically by comparing its prediction with the experimental or 
numerical-experimental NNS distributions for a number of mixed systems. We 
confine our consideration to the systems that conserve time reversal systems. 
Our purpose is to find out how far can we describe the Poisson-GOE transition 
in terms of Eq. ([20]) with (3—1. Explicitly, the K-deformed Wigner surmise for 
Poisson-GOE transition reads 

-Pk-goeK s) = A(k)s exp K [-B{k)s 2 ] , (29) 

where 

A(k) = 2(1 -k 2 )B(k), B(/s) = | 

As mentioned above, the requirement of finite second moment of the spacing 
distribution restricts the range of variation of k to be 

0<«<i (31) 

with k = yielding Pgoe(s)- The evolution of the distribution Pk_goe(1/2, s) 
in this interval is shown in Fig. 1. As k increases, the peak of -Pk-GOe(k, s) 
moves towards small value of s. At k = 1/2, the distribution has its peak at 
the same position as the semi-Poisson distribution 

P SP (s) = 4sexp(-2s). (32) 

The semi-Poisson distribution was suggested to describe a narrow intermediate 
region between insulating and conducting regimes exemplified by the Anderson 
localization model [33], with the two limiting cases being described by Poisson 
and Wigner statistics, respectively. It was introduced to mimic new seemingly 
universal properties in certain classes of systems, in particular, being character- 
istics of the "critical quantum chaos" . Figure 1 compares the limiting supersta- 
tistical distribution, Pk-goe(1/2, s), with the semi-Poisson distribution. We see 
from the figure that the two distributions are quite similar in shape. Therefore 
it is worth in the present context to use the semi-Poisson statistics as a reference 
distribution marking the limit of validity of the base-invariant random-matrix 
description of mixed systems. Analogous conclusion was reached in the analysis 
of the transition out of chaos using the superstatistical generalization of RMT 

In the following, we test how accurately the distribution Pk-goe(k, s) repro- 
duces the histogram data obtained in a number of numerical experiments. 



(i-- 2 )r(£- 



e B/a 



(30) 



5.1 Random points on fractals 

Sakhr and Nieminen [35] have recently noted that the NNS distribution of ran- 
dom points uniformly distributed on a line is given by the Poisson distribution 
l|28p . On the other hand he NNS distribution of random points uniformly dis- 
tributed on a plane is given by the Wigner surmise for GOE (j2"7| . Exploiting this 
observation, they introduce a model for a crossover transition between Poisson 
and Wigner statistics consisting of random points on a continuous family of self- 
similar curves with fractal dimensions between 1 and 2. This model is special 
since the intermediate statistics are described exactly by the Brody distribution 

m 

f*Brody(w, a) = a^s" exp {~b u s" +1 ) (33) 

with u) = d s — 1, which is frequently used in the phenomenological analysis 
of mixed systems. As a concrete example, they study point processes on the 
family of Koch fractals in 2-dimensional space. These fractals can be thought 
of as the attractors of a one-parameter family of iterated function systems. 
They show that NNS distributions of the random points undergo a continuous 
transition from Poisson to Wigner statistics as the similarity dimension of the 
fractals d s varies from 1 (line) to 2 (plane). Their results are shown in Fig. 2 
together with the best-fit curves calculated using K-deformed Wigner surmise 
of Eq. (|2U)) . Figure 3 shows the absolute deviation between the theoretical and 
numerical-experimental curves, which is defined by 

1 N 

A(«) = — ^ |-Pk-GOe(k, Si) - P CX pcrimont(Si)l ( 34 ) 
i=l 

The two figures show that the agreement is satisfactory for the nearly chaotic 
regime (the lower panels in Fig. 2), and gradually decrease as the distributions 
become less chaotic to become unsatisfactory as the system approaches integra- 
bility (the upper panels in Fig. 2). If we accept an absolute deviation below 
A = 0.05 as an indicator for satisfaction of the agreement with the data, which 
in our case corresponds to a fractional dimension (and thus Brody's parameter) 
nearly equal to d s — 1.46 (u> — 0.46). We note that the Brody distribution that 
best-fits the semi-Poisson distribution has a parameter oj = 0.41. 

The transition between integrability and chaos in this numerical experiment 
is monotonic. The degree of "chaoticity" is measured by the similarity dimension 
of the fractals d s and modelled in this paper by the Kaniadakis index k. In 
many systems, however, the degree of chaoticity changes in a complicated way 
as a system parameter is varied monotonically, and so the energy-level statistics 
will not undergo a direct transition from Poisson to Wigner. In the following 
subsection we consider a system which undergoes a Poisson- Wigner-Poisson 
transition. 

5.2 Elliptical stadium billiard 

Lopac et al. [37] investigated the dynamical properties of the elliptical stadium 
billiard, which is a generalization of the Bunimovich circular stadium billiard 
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[55] . In its dependence on two shape parameters 5 and 7 , this system reveals 
a rich interplay of intcgrable, mixed and fully chaotic behavior. Their results 
for the classical dynamics confirm the existence of a large fully chaotic region 
surrounding the straight line (5 = 1-7 that corresponds to the Bunimovich bil- 
liard. The histograms of NNS distributions reported in [37] are for the elliptic 
billiards with S = 7. For these billiards, the chaotic region of the shape param- 
eters is around 6 = 0.5. Fig. 4 shows the result of comparison between these 
histograms with the predictions of Eq. (|29p . The best-fit values of k together 
with the absolute deviations A are given in the figure. If we again accept the 
value of A = 0.05 to quantify the start of unsatisfactory agreement with the 
data, then we consider the K-deformed approach to be unsatisfactory outside 
the region of 0.4 < d s < 0.75. This shows the the K-deformed RMT provides a 
better description for mixed systems as the evolve approaching equilibrium or 
departing out of it. 

5.3 Boson systems with one- plus two-body random inter- 
action 

Random matrix ensembles (and in particular GOE) for m-particle system, in- 
volve interaction up to m-body in character and are dominated by the m-body 
interactions. In real many-particle systems such as atoms, nuclei, quantum dots 
etc. however, the inter-particle interaction is known to be only 1-body and 
2-body in character. Hence the concept of embedded ensemble of fc-body in- 
teraction (EGOE(fe)), was introduced [39]. Embedded ensembles are usually 
studied using only 2-body interactions, i.e., by allowing the single-particle en- 
ergies to be degenerate. Patel et al. [JU] have demonstrated that the EGOE(2) 
for dense interacting boson systems exhibit spectral fluctuations of GOE type. 
On the other hand, ensembles with only one-body interaction would show the 
Poisson character. Chavda et al. [41] presented numerical results for the study 
of the statistical properties of embedded ensembles of interacting Bosons when 
the interaction is taken as H = Hi + XH2 , so that the strength of the 2-body 
interaction can be varied by changing the value of A; H n here represents the 
n-body interaction. They examined the nature of the transition in spacing dis- 
tribution as the strength of the 2-body interaction is varied while keeping the 
nondegenerate but fixed single-particle energies They studied the nature of the 
spacing distribution as a function of strength A of the 2-body interaction. Their 
results are compared with with the predictions of Eq. (|2"9"]l in Fig. 5. We see 
again that In conclusion, we can see that the agreement between the histograms 
of the NNS distributions for bosonic ensembles and the K-deformed Wigner 
surmise steadily improves as the distribution moves from Poisson type to the 
Wigner type. 

5.4 Statistics of 2 + levels in even— even Nuclei 

Originally, RMT was proposed for the purpose of describing the spectral prop- 
erties of nuclear levels. Haq et al [22] have shown that the spectral fluctuation 



10 



properties of nuclear resonances near nucleon emission threshold are in very 
good agreement with the predictions of GOE. On the other hand, NNS dis- 
tributions of nuclear levels near the ground state lie between the Poisson and 
Wigner distributions [33]- Abul-Magd et al. [H] have collected the empirical 
information on low-lying levels of even-even nuclei with spin-parity 2 + , pub- 
lished in Nuclear Data Sheets until July 2002. To obtain statistically relevant 
samples, the nuclei are grouped into classes defined by the ratio R4/2 of the ex- 
citation energies of the first 4 + and 2 + levels. This ratio serves as a measure of 
collectivity in nuclei. Analyses done in Ref. [33] show that the NNS distribution 
for each group vary strongly with R4/2 and takes a more regular shape in nuclei 
that have one of the dynamical symmetries of the interacting Boson model [33] . 
Figure 6 shows the result of comparison between the histograms for each group 
with the predictions of Eq. ([29]) . The best- fit values of K together with the 
absolute deviations A are given in the corresponding windows. As seen in Fig. 
6, the entropic parameter k is indeed dependent on R4/2. It takes large value 
at i? 4 /2 = 2.0, 2.5, and 3.3, which respectively correspond to the U(5), SO(6), 
and SU(3) dynamical symmetries of the interacting Boson model. 



6 Transition intensity distribution for mixed sys- 
tems 

Matrix elements of transition operator probe the system's wave functions so that 
their statistical fluctuations provide additional information. In chaotic systems, 
the reduced transition probabilities follow the Porter- Thomas distribution [5J. 
The probability Bif of a transition from the initial configuration \i) to the final 
configuration |/) is given by 

B if = \W if \ 2 , (35) 

where Wif are the matrix elements of the transition operator in a special basis. 
If the transition operator conserves time reversibility, the matrix elements W% / 
are real. For a chaotic system with time reversal symmetry, it is reasonable 
to assume that Wif are identically-distributed Gaussian random variable. This 
entails that the transition intensities can be represented by a random variable 
that takes the values 

»" = (tf) <36) 

where (Bij) is a suitably defined local average value [46] . and has a Porter- 
Thomas distribution 

PMv) = J^<r m . (37) 
V n y 

This is a x 2 -distribution of one degree of freedom. The parameter 77 is defined by 
the requirement that (y) = 1, and is equal to 1/2. A more elaborate derivation 
of the Porter-Thomas distribution is given by Barbosa et al. [4"T] . 
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As the system becomes more regular, the transition probabilities deviate 
from the Porter-Thomas distribution. To account for these deviations, Alhassid 
and Novoselsky [35] suggested that the transition widths in mixed system may 
be analyzed in terms of a ^-distribution with v degrees of freedom 

p Ay^) = ¥ j^y v/2 -^- uv ' (38) 

When v = 1, This distribution coincides with the Porter- Thomas distribution. 
The /c-deformed distribution does not fit well the empirical distributions but 
consists better with the observed number of weak transitions as compared with 
the Porter-Thomas distribution (see, e.g. [3J3 133 133)- When the empirical 
distributions are expressed as functions of the log a y have peaks at log a y < 
while all the \ 2 distributions are peaked at log Q jy = for all values of a and 
v. We show in the present paper that the K-deformed RMT provides us with 
a generalization of the Porter- Thomas distribution, which is more suitable for 
the analysis empirical data. 



6.1 The /^-deformed Porter-Thomas distribution 

We now derive the K-deformed generalization of the Porter- Thomas distribu- 
tion. We do this by replacing the exponential function in (|3T[) by a K-deformed 
exponential. One then obtains 

P K -PT(y) = £ ^exp K [-d(«)y] (39) 

where the quantities c(k) and d(n) are obtained from the normalization condi- 
tion and the requirement that (y) = 1, which yield 

= 2 wjf_ + i!W ), (40) 

and 

p (-! — — ) r (— — \- — ) 

d(/s) - 4 ! 4 L - (41) 

We compare in Fig. 7 the evolution of Pk-ptO^v) during the transition from 
chaotic to regular dynamics by varying the tuning parameter k from (GOE) 
to 0.6 (almost regular) with a corresponding evolution of P x 2(\ny,v) where v 
varies between 1 and 0.1. Of special interest is the fact that the maximum of 
P % 2 (lny, v) occurs at lny = for any value of v as mentioned above. This prop- 
erty does not hold for Pk-pt (lny) . The peak of the ft-deformed Porter-Thomas 
distribution for less chaotic systems occurs at lny < and moves towards lower 
values as the parameter n increases. We show in the next subsection that this 
is indeed the behavior of physical systems. 
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6.2 Data analysis 

Hamoudi, Nazmitdinov and Alhassid [2H] calculated the electric quadrupole 
(E2) and magnetic dipole (Ml) transition intensities among the isospin T = 0, 1 
states of nuclei with mass number 60. They applied the interacting shell model 
with realistic interaction for pf shell nuclei with a 56 Ni core. It is found that 
the B(E2) transitions are well described by a GOE (Porter- Thomas distribu- 
tion). However, the statistics for the B(M1) transitions is sensitive to T z . The 
M2 transition operator consists of an isoscalar and isovector components. The 
T z = 1 nuclei, in which both components contribute, exhibit a Porter-Thomas 
distribution. In the meanwhile, a significant deviation from the GOE statis- 
tics for the T z = nuclei, where the matrix elements are purely isoscalar and 
relatively weak [50]. 

We use using the K-deformed Porter-Thomas distribution to analyze the 
reduced Ml transition intensities for both the T z = 1 60 Co nuclei and T z = 
60 Zn nuclei calculated by Hamoudi et al. H§] ■ These authors sampled a large 
number of matrix elements for each transition operator, which is equal to 56 2 = 
3136 and 66 2 = 4356. Figure 8 compares the results of calculations using Eq. 
(|3T)|) with the numerical results of Hamoudi et al. as well as the "best-fit" \ 2 
distribution deduced by these authors. The figure clearly shows the advantage of 
the using the K-deformed Porter-Thomas distribution over the x 2 distribution, 
at least for this numerical experiment. 

7 Conclusion 

Tsallis' entropy has been considered by several authors as a starting point for 
constructing a generalization of RMT. However, there are several other general- 
ized entropies. This paper derives K-deformed orthogonal and unitary random- 
matrix ensembles by maximizing the generalized entropy proposed by Kani- 
adakis under the constraints on normalization of the distribution and the ex- 
pectation value of Tr[H^H). We obtain a simple formula for the distribution 
function in the matrix-element space of both of the K-deformed GOE and GUE. 
This leads to new expressions for the NNS distributions of the eigenvalues when 
k ^ 0. The case of n = is the one that leads to the Wigner surmise when 
the entropy is given by the Shannon measure. The high accuracy of Wigner's 
distribution describing chaotic systems justifies the use of this distribution as a 
starting point for introducing the k deformations which follow as a consequence 
of adopting the Kaniadakix entropy. As in the case of Tsallis' entropy, the NNS 
distribution obtained here for Kaniadakis' entropy has an intermediate shape 
between the Wigner distribution and the Poissonian that describes generically 
integrable systems. However, neither of these distributions reaches the Poisso- 
nian form. We have explored the possibility to us the NNS distribution obtained 
in this paper for modelling systems with mixed regular-chaotic dynamics at least 
when they are not far from the state of chaos. We have tested these expres- 
sions by comparing their predictions with the NNS distributions obtained in 
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a number of numerical experiment. We also introduce a k deformation to the 
Porter-Thomas distribution that describes the transition intensities in chaotic 
systems. The resulting transition- intensity distribution is found to agree with 
the results of shell model calculation for a number of nuclei better that the % 2 
distribution, which is often used for this purpose. 
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Figure caption 



FIG. 1. NNS distributions obtained by using the Kaniadakis statistics. The 
entropic indices used are k = (GOE),0.3 and 0.5. The semi=Poissonian dis- 
tribution is also plotted to show that it has a peak nearly at the same position 
as that of the limiting case of the K-deformed Wigner surmise (with n = 1/2). 

FIG. 2. NNS distributions obtained in a point processes on the family of 
Koch fractals in a 2-dimensional space [35] are shown by histograms. The simi- 
larity dimension d s of the fractals are respectively d s = 1.002, 1.060, 1.140, 1.267, 1.375, 1.496, 1.684 , 1.812 
and 1.988. The continuos curves are calculated by Eq. ([29|i with values of the 
entropic index k shown in each window. 

FIG. 3. The fractal dimensional d s corresponding to the NNS distribu- 
tions shown in Fig. 2 plotted against the best-fit values of entropic index n. 
The continuos line is a parametrization of the empirical values by d s (n) — 
2.01-e~( K -°- 573)2 /( 2x0 - 1342 ). 

FIG. 4. NNS distributions of levels of the elliptical stadium billiards with 
shape parameters 7 = S [37], fitted to the distributions calculated by Eq. (|29l) 
with values of the entropic index k shown in each window. 

FIG. 5. NNS distributions for interacting bosons with Hamiltonian H = 
Hi + XH 2 , for various values of A calculated in [UJ (a) with 10 bosons in 4 
single-particle states and (b) with 10 bosons in 5 single-particle states, fitted 
to the distributions calculated by Eq. (I2TH) with values of the entropic index 
k shown in each window. The last histograms in both (a) and (b) are for 
EGOE(2). 

FIG. 6 NNS of 2 + of even-even nuclei in different ranges of the R4/2 ratio 
fitted to the distributions calculated by Eq. (f29| with values of the entropic 
index k shown in each window. 

FIG. 7. Evolution of the K-deformed Porter-Thomas distribution f^.px(lny) 
and the \ 2 distribution P x 2(lnj/, v) during the transition from chaotic to regular 
dynamics. The solid curves, labeled as PT, refer to the Porter-Thomas distri- 
bution. 

FIG. 8. Nuclear shell-model Ml transition intensities in A = 60, calculated 
by Hamoudi et al. |49] . (histograms) compared with the K-deformed Porter- 
Thomas distribution (|3"9"|) with parameters k = 0.38, 0.61 and 0.63, respectively 
(solid curves) and the \ 2 distribution (17) with parameters v = 1, 0.64 and 0.34, 
respectively (dashed curves). 
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